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PREFACE 

This  report  was  prepared  by  Dr.  Donald  E.  Nevel,  Research  Physical  Scientist,  Applied  Research 
Branch,  Experiniental  Engineering  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory.  The  report  was  published  under  DA  Task  1T0621 12A13001,  Cold  Regions  Research  — 
Applied  Research  and  Engineering. 

The  author  wishes  to  thank  Ward  MacKenzie,  who  made  the  numerical  evaluation,  including 
writing  of  the  detailed  computer  program. 
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CONVERSION  FACTORS:  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 

These  conversion  factors  include  all  the  significant  digits  given  in 
the  conversion  tables  in  the  ASTM  Metric  Practice  Guide  (E  380), 
which  has  been  approved  for  use  by  the  Department  of  Defense. 
Converted  values  should  be  rounded  to  have  the  same  precision  as 
the  original  (see  E 380). 


Muitipiy  By  To  obtain 

inch  25.4*  millimeter 


BEARING  CAPACITY  OF  RIVER  ICE  FOR  VEHICLES 


Donald  E.  Nevel 


INTRODUCTION 

In  cold  environments  vehicles  may  cross  a river  over  the  ice.  If  the  river  ice  is  not  thick  enough 
to  support  the  vehicles,  an  ice  bridge  may  be  built  by  flooding  the  surface  of  the  ice  with  a few 
inches  of  water  and  waiting  for  the  water  to  freeze.  By  repeating  this  procedure,  the  thickness  of 
the  ice  can  be  increased  faster  than  by  the  natural  growth  of  ice.  This  procedure  should  not  be  re- 
peated indefinitely  since  there  is  a maximum  ice  thickness  which  can  be  obtained.  This  maximum 
thickness  depends  upon  the  local  climatic  environment.  If  snow  is  on  the  ice  it  should  not  be  re- 
moved. However,  if  the  snow  depth  is  greater  than  5 in.,  the  snow  should  be  compacted  before 
being  flooded. 

In  the  past,  two  types  of  ice  bridges  have  been  built.  For  one  type,  the  ice  was  flooded  within 
a confined  strip  from  shore  to  shore.  For  the  other,  the  ice  was  flooded  freely  from  shore  to  shore. 
G.  Frankenstein  of  CRREL  has  shown  that  free  flooding  is  easier  to  do  and  produces  a stronger 
bridge  than  confined  flooding.  This  was  demonstrated  during  a school  for  crossing  river  ice  given 
at  Fort  Wainwright,  Alaska,  in  November  1965  by  G.  Frankenstein,  R.  Garner  (also  of  CRREL), 
and  the  author.  Hence,  it  is  expected  that  future  ice  bridges  will  be  built  by  free  flooding. 

In  a previous  report  by  the  author  (Nevel  1965b),  equations  describing  the  stresses  in  an  i(e 
bridge  built  by  confined  flooding  were  presented.  Although  these  equations  will  probably  not  be 
used  for  practical  applications,  eq  1 5c  and  1 5d  of  this  1965  report  are  incorrect.  These  equations, 
respectively,  should  be 


0^+a^  =5^+0^ 

dy^  dx^  y=j  dy^  Bx‘ 

O ^+(2 -a)  = D ^ 

by^  3x2  dy  ay  3 


3x2 


0 21Sil  (2  - o) 

3y^  3x2  ay 


+ (2-a) 


3^^  ay  y=0 


Additional  equations  which  result  from  these  two  equations  are  also  incorrect. 

For  the  bridge  built  by  free  flooding,  the  ice  thickness  varies  but  its  variations  are  not  precisely 
known.  Free  flooding,  however,  produces  a very  wide  bridge  which  may  be  approximated  by  an 
infinite  width  bridge.  Therefore,  for  practical  applications  the  bearing  capacity  of  a river  that  has 
a uniform  ice  thickness  will  be  determined.  The  result  will  then  represent  the  bearing  capacity  of 
the  frozen  river  without  an  ice  bridge  or  with  an  ice  bridge  built  by  free  flooding.  Because  of  the 
formation  of  shoreline  cracks,  a simple  support  is  the  best  representation  of  the  boundary  conditions 
along  each  shore  of  a river.  Hence,  we  want  to  analyze  an  infinite  strip  on  an  elastic  foundation 
simply  supported  on  opposite  sides. 


THEORY 


} 


Solution  by  superposition 

A.D.  Kerr  (1959)  used  the  method  of  images  to  solve  a number  of  floating  ice  sheet  problems 
which  have  simply  supported  boundary  conditions.  One  of  his  solutions  considers  a simply  sup- 
ported strip  on  an  elastic  foundation  with  a concentrated  load.  His  method  is  easily  extended  for 
distributed  loads. 

Fur  a single  distributed  load  on  an  infinite  ice  sheet,  the  deflection  for  a point  not  under  the  load 
is  given  by  M.  Wyman  (1950)  as 

w{R)  = -£-  i |ber'(A/8)  ker«-bei'(b/K)  kei/?l  (1) 


where  P 
ft 
b 

8^ 

R 

r 

h 

E 

a 


total  load 
density  of  water 
radius  of  loading 

-o^) 

r/8 

disunce  from  the  center  of  loading  to  the  point  under  discussion 
ice  thickness 
modulus  of  elasticity 
Poisson's  ratio. 


Q>nsider  an  infinite  ice  sheet  with  two  distributed  loads,  one  acting  downward  and  one  acting 
upward.  The  perpendicular  bisector  of  the  line  joining  these  two  loads  is  a line  of  zero  deflection. 
The  bending  moments  along  this  line  are  also  zero.  Hence  this  perpendicular  bisector  acts  as  a 
simple  support  in  the  ice  sheet.  This  is  demonstrated  in  Figure  1 where  one  load  acts  downward 
(positive)  at  Xq,  and  another  load  acts  upward  (negative)  at  -Xq,  y^.  The  perpendicular  bisector 
is  they  axis  which  acts  as  a simple  support. 

In  order  to  make  the  linex  = Z.  a simple  support  in  Figure  1,  an  infinite  set  of  these  positive 
and  negative  loads  must  be  placed  on  the  ice  sheet.  The  first  two  such  sets  are  shown  in  Figure  1 . 

Hence  by  superposition,  the  deflection  at  point  Q is  given  by 


M^,y)  = 53  IM^n)  - 


(2) 


where 

R„^  = (x-Xo-2Z.n)2  + (v-yo)2 

= (>f  + JtO-2^n)^  + (y->'o)^- 
The  angles  associated  with  these  distances  are 


tan 


0=_ZZ>:o 


x-Xq-2L„ 


x+Xo-2Z.„ 


The  bending  moments  can  be  easily  determined  by 

0(1  g)  / d^w  . J_  \ 

2 RdR/ 


2 


(3) 


(4) 


2 


= -M. 


■ g)  /d^w  J_  dw  \ 

la/?2  " « a/?/ 


W = - 0(1  - °)  _ 1 |^)5in20 

2 \ 3/;2  R dR  / 


cos  20 


(5) 

(6) 


where  R is  any  of  the  R„  or  R„  and  the  0 the  corresponding  6 or  S. 

After  the  appropriate  sum  is  obtained,  the  maximum  moment  may  be  obtained  by  using  Mohr’s 
circle.  The  maximum  stress  is  then 

If  the  point  is  under  the  load  (r  < b),  the  corresponding  term  in  the  series  should  be  replaced 
by 


w(R)  = — + |.  [ker'(b/8)  ber/?  - kei'(b/6)bei/?]l  (7) 

7r/f82  I /,2  ft  I 

However,  eq  7 is  not  valid  for  determining  moments  if  the  load  is  highly  concentrated.  A three 
dimensional  theory  of  elasticity  has  been  given  (Nevel  1970)  and  evaluated  forr  = 0.  Numerically 
these  results  agree  with  a formula  proposed  by  Westergaard  (1926).  When  b/Zt  > 1 .724,  Westergaard 
says  that  cq  8,  which  is  a direct  result  of  eq  7,  should  be  used. 


= /My  = Ij:  (1 +a)|  kei'(b/8) 

^xy  = 0. 

When  b/b  < 1 .724,  b in  eq  8 should  be  replaced  by  a where  o = (1 .6b2  + b2)’'4  _ 5b/8. 


(8) 


Solution  by  Fourier  integral 

The  solution  for  a load  uniformly  distributed  over  a rectangular  area  on  a river  has  been  obtained 
(Novel  1965a).  The  solution  of  the  differential  equation  was  obtained  by  taking  a Fourier  integral 
transform  along  the  river  followed  by  a Laplace  transform  across  the  river.  The  inverse  Laplace 


transform  was  obtained  but  the  inverse  Fourier  integral  transform  could  not  be  obtained.  Thus  the 
solution  is  expressed  in  integral  form.  This  solution  will  not  be  used  because  in  the  following  section 
a simpler  solution  is  obtained  for  the  same  problem.  However,  the  advantage  of  the  Fourier  integral 
solution  is  that  boundary  conditions  other  than  simple  supports  may  be  met  at  each  shore. 


Solution  by  Fourier  series 

The  solution  for  a uniform  load  distributed  over  a rectangular  area  on  an  ice  bridge  has  been  ob- 
tained (Nevel  1965b).  The  solution  of  the  differential  equation  was  obtained  by  taking  a finite 
Fourier  transform  across  the  river  followed  by  a Laplace  transform  along  the  river.  The  inverse 

Laplace  transform  was  obtained;  this  left  the 


RIVER  ICE 


Figure  2.  Concentrated  loads. 


solution  in  series  form.  By  requiring  this  solu- 
tion to  meet  different  boundary  conditions, 
the  solution  of  the  river  ice  problem  is  ob- 
tained. 

Thus  for  two  concentrated  loads  on  the 
ice  as  shown  in  Figure  2,  we  require  that  all 
of  equations  8 in  Nevel  (196Sb)  go  to  zero 
as  T)  approaches  infinity.  This  gives  A = 
(e'^/e2)(|3s  + 7c)„  and  8 = e'<’“(^ls  - yc)„ 
(see  Nevel  (1965b)  for  all  notations] . Thus 
the  solution  is 


= -Ar  f ^ /sin^Tsin^f 


Afy  = -/>  J £ (y"'-a02  V')sin0rsin«? 


on 

(a/"  - 0^  y)  sin 01’ sin 05 


^xy  = -^(1  -ff)  r £ 0'''5i'’0rcos0{ 


y = (13s  + 7c)  -t-  W(o  - 17)  IPs  + 7c)„.^  + 

2e^  2€'‘ 


+ W(»?  - a)  ^ (/3s  x)„-a 

2e'' 


y = —^—2 5a-„-W(T}-a)-S^-5 


+ H(t,  - a)  ‘ IPs  - Trln-o 
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Figure  3.  Distributed  loads. 


For  two  loads  distributed  over  a rectangular  area  as  shown  in  Figure  3,  we  require  that  all  eq  10 
in  Nevel  (1965b)  go  to  zero  as  rj  approaches  infinity.  This  gives 


e*  I I 


Thus  the  solution  is 


w = T S sin0rsin(^| 

^ ^ ‘t»'F  ^ 


M,.  = - 


^ f sin-^  sin0rsin0| 


M, 


= {oy^  -<I>^  y)  sin.^  sin 0r  sin 01 


t=l 


/Mxy  = -/’(1-ff)|  2 sin-^  sin0rcos0J 


t=i 


where 


/,  = - -L- 

L-ax  _ j_ 

jw(x)(e-P''c-1) 

a+#</2-T) 

2e'‘ 

1 J V*a-ul2  2e* 

a-nl2-n 

+ J_ 

j/y(x)(e-^c-1)j 

r\-a-nl2 

2e'‘ 

T)-a+(i/2 

■'i " ^ -i- 


a+M/2-r) 

■M/2-rj 


(13 


(14 


(15 


(16 


y 


n 

1 


1 

2 


r)+o+ji/2 

T7+a-M/2 


1 

2 


H(x)  e-^  s 


a+»i/2-i7 

a-ul2-ri 


+ 


+ 


1 

2 


tj-<»-»i/2 

Ti-a+M/2 


Here  again  the  above  moments  are  not  valid  for  relatively  concentrated  loads  if  the  point  under 
discussion  is  under  .he  load.  A three  dimensional  theory  of  elasticity  solution  is  not  readily  avail- 
able for  a recungular  load.  Therefore,  if  the  load  is  relatively  concentrated,  the  rectangular  load 
should  be  replaced  wtih  an  equivalent  circular  load  and  Westergaard’s  formula  should  be  used. 


APPLICATION  AND  RESULTS 

The  bearing  capacity  of  river  ice  for  the  following  vehicles  is  given  in  Appendix  A,  Figures  Al- 
All. 

Tracked  vehicles 

1.  M48 

2.  M8 

3.  SP  M109  How 

4.  PCM113 

5.  PCM577 

Wheeled  vehicles 

6.  M51,  5-ton  dump  truck  w/winch  and  3J5-ton  Bolster  trailer 

7.  M52,  5-ton  tractor  w/winch  and  Ml  72A1 , 25-ton  low  bed  semi-trailer 

8.  20-ton  crane. 

The  detailed  specifications  of  these  vehicles  are  given  in  Appendix  B. 

For  the  tracked  vehicles,  it  is  obvious  that  the  solution  by  Fourier  series  should  be  used  since  the 
contact  area  of  the  track  is  rectangular.  Each  vehicle  was  placed  at  the  center  of  the  river  since  this 
appears  to  be  the  most  critical  position.  The  moments  were  obtained  under  the  center  of  the  track 
and  then  the  maximum  stress  was  obtained.  The  results  are  shown  in  Figures  A1-A4. 

The  M51  and  M52  wheeled  vehicles  have  12-ply  1 1 x20  tires  inflated  to  70  psi.  For  these  tires 
Freitag  and  Green  (1962)  have  shown  that  the  contact  area  more  nearly  represents  a rectangle  than 
a circle.  Hence,  the  Fourier  series  solution  will  be  used  for  these  vehicles.  The  Fourier  scries  solu- 
tion will  also  be  used  for  the  remaining  wheeled  vehicles  for  convenience  since  no  contact  area  was 
available  for  the  tires  of  these  vehicles.  The  contact  area  of  each  tire  was  obtained  by  dividing  the 
load  on  the  tire  by  the  inflation  pressure.  The  ratio  of  length  to  width  of  the  contact  area  of  tires 
varies  from  about  1.5  to  2.  The  ratio  of  1.5  was  selected  since  this  should  produce  higher  stresses. 

The  stresses  at  the  center  of  the  contact  area  of  each  tire  on  each  vehicle  were  determined  as 
follows.  The  tire  was  placed  at  the  center  of  the  river  since  this  seems  to  be  the  most  critical  posi- 
tion. Westergaard’s  formula  was  used  to  evaluate  the  moments  due  to  the  load  of  this  tire  since  the 
contact  area  is  relatively  concentrated.  The  additional  moments  due  to  the  remaining  tire  loads 
were  evaluated  by  the  Fourier  series  solution  and  added  to  the  moments  previously  obtained  from 
Westergaard’s  formula.  The  maximum  stress  under  the  tire  being  tested  was  then  obtained  by  using 
Mohr’s  circle.  A comparison  of  the  maximum  stress  under  each  tire  was  made  to  determine  the 
critical  stress  of  the  vehicle.  The  results  are  shown  in  Figures  A5-A1 1 . The  Fortran  II  computer 
program  for  the  computations  is  given  in  Appendix  C. 

For  the  vehicles  considered,  the  stress  in  kilograms  per  square  centimeter  divided  by  the  total 
load  in  tons  (2000  lb)  is  presented  in  Figures  A1-A1 1 as  a function  of  the  width  of  the  river  and 
the  ice  thickness.  The  total  load  means  the  weight  of  the  vehicle  plus  the  load  that  is  being  carried. 
For  trucks,  the  graphs  for  completely  loaded  and  unloaded  conditions  arc  presented.  For  a lightly 
loaded  (up  to  25%)  truck,  the  graph  for  the  unloaded  truck  should  be  used.  For  any  other  partially 
loaded  (over  25%)  truck  the  graph  for  the  loaded  truck  should  be  used.  The  allowable  stress  for 
cold,  clear,  fresh  water  ice  (as  can  be  expected  in  Alaska)  is  10  kg/cm’.  If  the  ice  is  warm  or  partially 


6 


deteriorated,  the  allowable  stress  is  S kg/cm’.  In  the  spring  thaw,  the  allowable  stress  can  drop  as 
low  as  1 to  2 kg/cm’ . 


As  an  example,  in  using  these  graphs,  suppose  we  have  a 75-ft-wide  river  covered  with  20-in.- 
thick  ice.  Wc  desire  to  take  across  the  ice  an  M51  truck  that  is  half  loaded.  Referring  to  Figures 
A5  and  A6,  the  truck  weighs  1 1.3  tons  and  has  a carrying  capacity  of  10.7  tons.  This  means  that 
the  total  weight  of  the  half  loaded  truck  is  16.7  tons.  We  select  Figure  A6  because  the  truck  is 
more  than  25%  loaded.  For  a 75-ft  width  and  20-in.-thick  ice,  we  find  that  Stress/Total  Load 
= 0.44.  The  stress  in  the  ice  will  be  0.44x(rof<7/ /.ootf)  = 0.44x16.7  = 7.5  kg/cm*.  This  is  a 
permissible  value  of  stress  for  cold,  clear,  fresh  water  ice,  but  not  for  partially  deteriorated  ice. 

In  these  graphs,  the  stress  factor  depends  upon  the  ice  thickness  and  the  width  of  the  river.  In 
order  to  eliminate  the  width  of  the  river,  we  choose  the  width  of  the  river  that  gives  the  maximum 
stress  factor.  For  these  critical  river  widths,  the  stress  factor  depends  upon  the  ice  thickness  as 
shown  in  Figure  A12.  If  the  river  is  very  wide,  the  stress  factor  does  not  depend  upon  the  river 
width  but  on  the  ice  thickness  as  shown  in  Figure  A13. 

In  summary.  Figure  A1 2 can  be  used  as  the  basis  for  the  bearing  capacity  of  river  ice.  For  a 
more  exact  analysis  Figures  A1-A1 1 can  be  used.  A speed  well  below  the  critical  velocity  should 
be  maintained. 
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WIDTH  OF  RIVER  IN  FEET 


ICE  THICKNESS  IN  INCHES 


Total  weight 
(tons) 


1) 

PCM  113  or  PCM577 

11.5 

2) 

M8 

21.1 

M8  with  commo  van  (AN/MR069) 

28.2 

3) 

20>ton  crane 

28.5 

SPM109  how 

26.2 

4) 

M51  5-ton  dump  truck  with  winch  - loaded 

22.0 

5) 

M51  5-ton  dump  truck  with  winch  - unloaded 

11.3 

6) 

M51  5-ton  dump  truck  with  winch  - loaded 
and  3.5-ton  bolster  trailer  - loaded 

26.2 

7) 

M52  5- ton  tractor  with  winch 

9.5 

8) 

M52  5-ton  tractor  with  winch  and 

M172  A1  25-ton  semi-trailer  - unloaded 

16.9 

9) 

M52  5-ton  tractor  with  winch  and 

M172  A1  25-ton  semi-trailer  - loaded 

42.5 

10) 

M48  tank 

49.0 

Figure  A 12. 


20  24  28 

ICE  THICKNESS  IN  INCHES 


Total  weight 
(tons) 


1)  PCM113or  PCM577 

2)  MS 

MS  with  coramo  van  (AN/MRC-69) 

SPM109  how 
20'ton  crane 

3)  M51  5-ton  dump  truck  with  winch  - loaded 

4)  M51  5-ton  dupip  truck  with  winch  - unloaded 

5)  M52  5-ton  tractor  with  winch 

6)  M5 1 5-ton  dump  truck  with  winoh  •*  loaded 
and  3.5-ton  bolster  trailer  - loaded 

7)  M52  5-ton  tractor  with  winch  and 

M172  A1  25- ton  low  bed  semi- trailer  - unloaded 
S)  M52  5-ton  tractor  with  winch  and 

M172  A1  25-ton  low  bed  semi  trailer  - unloaded 
9)  M4Stank 
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